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On the basis of frequentist analyses of experimental constraints from electroweak precision data, {g — 2)^, B 
physics and cosmological data, we investigate the parameters of the constrained MSSM (CMSSM) with universal 
soft supersymmetry-breaking mass parameters, and a model with common non-universal Higgs masses (NUHMl). 
We present likelihood functions for the masses of supersymmetric particles and Higgs bosons, as well as 
BR(6 sj), BR(_Bs fi~) and the spin-independent dark matter scattering cross section, o-p'. In the 

CMSSM we find preferences for sparticle masses that are relatively light. In the NUHMl the best-fit values for 
many sparticle masses are even slightly smaller, but with greater uncertainties. The likelihood functions for most 
sparticle masses are cut off sharply at small masses, in particular by the LEP Higgs mass constraint. Both in 
the CMSSM and the NUHMl, the coannihilation region is favoured over the focus-point region at about the 3-a 
level, largely but not exclusively because of {g — 2),^. Many sparticle masses are highly correlated in both the 
CMSSM and NUHMl, and most of the regions preferred at the 95% C.L. are accessible to early LHC running, 
though high-luminosity running would be needed to cover the regions allowed at the 3-a levels. Some slepton 
and chargino/neutralino masses should be in reach at the ILC. The masses of the heavier Higgs bosons should be 
accessible at the LHC and the ILC in portions of the preferred regions in the (Ma, tan /3) plane. In the CMSSM, 
the likelihood function for BR(_Bs fi^ jj,^) is peaked close to the Standard Model value, but much larger values 
are possible in the NUHMl. We find that values of > 10~^° pb are preferred in both the CMSSM and the 
NUHMl. We study the efi'ects of dropping the (g — 2)^, BR(6 sj), Q.^h? and Mh constraints, demonstrating 
that they are not in tension with the other constraints. 
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1. Introduction 

Supersymmetry (SUSY) [1~3] is one of the 
favoured ideas for physics beyond the Stan- 
dard Model (SM) that may soon be explored at 
the Large Hadron Collider (LHC). In a recent 
paper [4], we presented some results from fre- 
quentist analyses of the parameter spaces of the 
constrained minimal supcrsymmetric extension of 
the Standard Model (CMSSM) — in which the 
soft supersymmetry-breaking scalar and gaugino 
masses are each constrained to universal values 
mo and mi/ 2, respectively [5-24] — and the 
NUHMl — in which the soft supersymmetry- 
breaking contributions to the Higgs masses are 
allowed a different but common value [25-27]. 
Other statistical analyses in these models can be 
found in [28-39] and Markov Chain Monte Carlo 
(MCMC) analyses in [40-57]. For comparison, 
see also [58,59] for recent analyses in the next-to- 
minimal extension of the SM, as well as [60-62] for 
other analyses in supcrsymmetric models without 
a dedicated fit. 

The results presented in [4] included the param- 
eters of the best-fit points in the CMSSM and the 
NUHMl, as well as the 68 and 95% C.L. regions 
found with default implementations of the phc- 
nomcnological, experimental and cosmological 
constraints. These include precision electroweak 
data, the anomalous magnetic moment of the 
muon, [g — 2)^, _B-physics observables (the rates 
for BR(6 S7) and BR(_B„ tv^), Bg mix- 
ing, and the upper limit on BR(i?s ^ 
the bound on the lightest MSSM Higgs boson 
mass, Mh, and the cold dark matter (CDM) den- 
sity inferred from astrophysical and cosmological 
data ^, assuming that this is dominated by the 
relic density of the lightest neutralino, Vl^h? . We 
also discussed in [4] the sensitivities of the areas 
of the preferred regions to changes in the ways 
in which the major constraints arc implemented. 
We found that the smallest sensitivity was to the 
CDM density, and the greatest sensitivity was 
that to {g — 2)^. 



^We did not include the constraint imposed by the exper- 
imental upper limit on the spin-independent DM scatter- 
ing cross section cr^', which is subject to astrophysical and 
hadronic uncertainties, as discussed below. 



In this paper we adopt the frequentist approach 
from [4], which is different from the Bayesian 
approach adopted in [41-52]. A key issue in a 
Bayesian approach is the appropriate choice of 
priors. As discussed in some recent Bayesian 
analyses of the CMSSM [45,46,52], conclusions 
for preferred regions of parameter space can de- 
pend the choice of priors. In our view, the re- 
sults of a Bayesian approach should not be con- 
sidered definitive unless they are shown to be suf- 
ficiently independent of plausible variations in the 
choice of priors. In our frequentist analysis, we 
use the MCMC technique to sample efficiently the 
CMSSM and NUHMl parameter spaces, and we 
generate sufficiently many chains to sample ad- 
equately these parameter spaces, as discussed in 
more detail in Section 2 of this paper. 

Our treatments of the experimental constraints 
from electroweak precision observables, _B-physics 
observables and cosmological data are, in general, 
very similar to those in [4]. Accordingly, we do 
not discuss details in this paper, contenting our- 
selves with a brief recapitulation and update. 

In Section 3 we extend the presentation of re- 
sults from our MCMC frequentist analysis to in- 
clude the global likelihood functions for var- 
ious observables, including Mh, BR(6 — > S7), 
BR(i3s fi~^fj,~) and the spin-independent DM 
scattering cross section, tTp\ as well as sparti- 
cle masses. We also discuss the correlations be- 
tween pairs of these observables, and compare the 
results in the CMSSM and NUHMl. We pay 
particular attention to the prospects for detect- 
ing SUSY in forthcoming experiments, including 
searches at the LHC and the ILC as well as B 
physics and direct searches for CDM. 

We present an update on the prediction of Mh 
in the CMSSM [53] and the first prediction for 
Mh in the NUHMl. For these analyses the ex- 
perimental constraints on Mh itself have been left 
out of the fit. The result in the CMSSM of [53] 
is confirmed with a best-fit value slightly below 
the LEP bound. Within the NUHMl, however, a 
value above the LEP bounds arises naturally. For 
other observables, however, the Mh information 
is included in calculating the likelihood func- 
tions. The likelihood functions for generic sparti- 
cle masses are skewed, being cut off at low masses 
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by the LEP lower limit on My^, in particular. On 
the other hand the likelihood functions rise more 
gradually for large masses, with the largest contri- 
bution arising from (g — 2)^. We see that the role 
of the Mh constraint is smaller in the NUHMl 
than in the CMSSM, reflecting the fact that the 
other constraints suggest, in the NUHMl, a value 
oi Mh somewhat larger than the LEP lower limit. 

As remarked in [4] , the preferred values of the 
sparticle masses are generally somewhat lower in 
the NUHMl than in the CMSSM. This is be- 
cause the extra degree of freedom in the Higgs 
sector allows lower values of 7711/2 to be reconciled 
with upper limits on deviations from the SM and 
the LEP lower limit on Mh- Recall that in the 
CMSSM, the Higgs mass mixing parameter, /i, 
and the Higgs pseudoscalar mass, Ma, are fixed 
by the minimization of the Higgs potential ensur- 
ing electroweak symmetry breaking when tan (3 
is chosen as an input parameter. In contrast, 
in the NUHMl, either ^ or Ma can be chosen 
as an additional input parameter ^, thus allow- 
ing substantial additional freedom in the light 
Higgs scalar mass for a given set of CMSSM pa- 
rameters (mi/2, ^0, ^0, tan/3). The greater free- 
dom in the Higgs sector also results in different 
mass ranges being favoured for the heavier Higgs 
bosons A, and for the heavier neutralinos, 
as observed in [4] . 

We find here that sparticle masses are mostly 
highly correlated. This could be expected for 
rrij^o and m^, which are both determined essen- 
tially uniquely by mi/2- However, the correla- 
tion is only slightly weakened for the slepton and 
squark masses, including m^^ . This is partly be- 
cause the largest contributions to the preferred 
values of most of these particles arc due to mi/2, 
rather than to mg. This tendency is reinforced 
by the fact that our likelihood analysis finds that 
the coannihilation regions are favoured in both 
the CMSSM and the NUHMl. However, this 
preference is slightly weakened in the NUHMl, 
where direct-channel annihilation through the 
heavy Higgs {A,H) poles may also play a sub- 
sidiary role, and larger values of toq become pos- 

■^The choice of either or Ma as an input is equivalent to 
a choice of the soft Higgs mass rrif^^ = ruf^^ mo at the 
GUT scale. 



sible. In particular, the correlation between m^j 
and TTig is particularly weak in the NUHMl, re- 
flecting the appearance of preferred regions of the 
parameters away from the coannihilation strip. In 
general, there are good prospects for discovering 
SUSY in early LHC running, in both the CMSSM 
and the NUHMl. 

We find that BK{Bs fi^fj,^) is expected to 
be close to its SM value in the CMSSM, because 
of the strong preference for relatively low tan/3 
where the supersymmetric contributions to this 
channel are small. They may be much larger in 
the NUHMl because of the freedom to choose 
Ma below its nominal CMSSM value. Spin- 
independent scattering of supersymmetric dark 
matter may well be observable in planned ex- 
periments in both the CMSSM and the NUHMl, 
where a somewhat larger range for is preferred 
in the NUHMl [26,63]. 

However, these optimistic conclusions rely crit- 
ically on the implementation of the [g — 2)^ con- 
straint using low-energy e+e^ data, as used in 
our analysis, and we discuss in Section 6.1 the 
implications of removing the {g — 2)^ constraint. 
We also discuss the predictions of our fits for 
BR(6 57), ri^/i^ and Alh, presenting the likeli- 
hood functions for each of these observables with- 
out their own contributions. None of these ob- 
servables exhibits any significant tension with the 
others. 

2. Description of the Frequentist Statisti- 
cal Method Employed 

We define a global likelihood function, which 
combines all theoretical predictions with experi- 
mental constraints: 

2 ^ (Ci - Pj)'^ 

+ x^(SUSY search limits) 

M I j-obs _ f fit 'i2 

Here N is the number of observables studied, 
Ci represents an experimentally measured value 
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(constraint) and each Pi defines a prediction for 
the corresponding constraint that depends on the 
supersymmetric parameters. The experimental 
uncertainty, a{Ci), of each measurement is taken 
to be both statistically and systematically inde- 
pendent of the corresponding theoretical uncer- 
tainty, a{Pi), in its prediction. We denote by 
X^iMh) and x^(BR(Bs n^i)) the contri- 
butions from the two measurements for which 
only one-sided bounds are available so far, as dis- 
cussed below. Furthermore we include the lower 
limits from the direct searches for SUSY parti- 
cles at LEP [64] as one-sided limits, denoted by 
"X^(SUSY search limits)" in eq. (1). 

We stress that, as in [4,53], the three stan- 
dard model parameters /sm = {^cthndi'mt, Mz} 
are included as fit parameters and allowed to 
vary with their current experimental resolutions 
^(/sm)- Wc do not include as a fit parame- 
ter, which would have only a minor impact on the 
analysis. 

Formulating the fit in this fashion has the 
advantage that the probability, P(x^,7Vdof), 
properly accounts for the number of degrees 
of freedom, A^dof, in the fit and thus repre- 
sents a quantitative and meaningful measure for 
the "goodness-of-fit." In previous studies [53], 
P(x^,A^dof) has been verified to have a flat dis- 
tribution, thus yielding a reliable estimate of the 
confidence level for any particular point in pa- 
rameter space. Further, an important aspect of 
the formulation is that all model parameters are 
varied simultaneously in the MCMC sampling, 
and care is exercised to fully explore the multi- 
dimensional space, including possible interdepen- 
dencies between parameters. All confidence levels 
for selected model parameters are performed by 
scanning over the desired parameters while min- 
imizing the function with respect to all other 
model parameters. That is, in order to determine 
the function x^(a;) for some model parameter x, 
all the remaining free parameters are set to values 
corresponding to a new minimum determined 
for fixed x. The function values where x^i^) is 
found to be equal to Xmin + determine the 
confidence level contour. For two-dimensional pa- 
rameter scans we use Ax^ ~ 2.28(5.99) to deter- 
mine the 68%(95%) confidence level contours. 



Only experimental constraints are imposed 
when deriving confidence level contours, with- 
out any arbitrary or direct constraints placed on 
model parameters themselves."^ This leads to ro- 
bust and statistically meaningful estimates of the 
total 68% and 95% confidence levels, which may 
be composed of multiple separated contours. Fi- 
nally, the sensitivity of the global fit to different 
constraint scenarios can be studied by removing 
one of the experimental constraints or by rescal- 
ing one of the experimental uncertainties, as dis- 
cussed in Sect. 3 in [4]. Studies of such scenar- 
ios are particularly helpful in identifying which 
experimental data are most useful in constrain- 
ing the theoretical model and hence in precisely 
studying how hyper-volumes in parameter space 
become more tightly constrained (either now or 
in the future). 

Since each new scenario in which a parameter 
is removed or an uncertainty re-scaled represents, 
fundamentally, a new function which must 
be minimized, multiple re-samplings of the full 
multi-dimensional parameter space are, in princi- 
ple, required to determine the most probable fit 
regions for each scenario. However, these would 
be computationally too expensive. To avoid this 
difficulty, we exploit the fact that independent 
X^ functions are additive and result in a well de- 
fined x^ probability. Hence, we define "loose" x^ 
functions, xfoosc which the term representing 
some constraint, e.g., r^cDM, is removed from the 
global x^ function. The xLoso function represents 
the likelihood that a particular set of model pa- 
rameter values is compatible with a sub-set of the 
experimental data constraints, without any ex- 
perimental knowledge of the removed constraint. 

An exhaustive, and computationally expen- 
sive, 25 million point pre-sampling of the Xkiose 
function is then performed in the full multi- 
dimensional model parameter space using a 
MCMC. Constraint terms representing the var- 
ious experimental scenarios are then re-instated 



^For reasons of stability of higher-order contributions, we 
limit the range of tan /3 to values below tan /3 = 60. As 
explained in Section 3 below, we furthermore impose a 
cut on parameter regions where the higher-order correc- 
tions relating the running mass to the on-shcU mass of the 
pseudo-scalar Higgs boson get unacccptably large. 
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or removed to form different functions, one 
for each scenario studied. If the scenario requires 
an additional constraint to be removed from the 
xfoosc function, the density of points pre-sampled 
for the x^oosg function was carefully tested and 
verified to also be an unbiased and sufficiently 
complete sampling of the studied model parame- 
ter space for the full function by using dedi- 
cated MCMC samples of approximately one mil- 
lion sampling points each, where the particular 
constraint in question was removed. Specifically, 
we use this technique to study the effects of re- 
moving individually the (g — 2)^, BR(6 — > 57), 
fi^ft,^ and Mh constraints. The precise values of 
the most probable fit parameters are determined 
via a full MINUIT minimization of the for each 
different scenario, but are performed only within 
the general parameter space regions not already 
excluded from the pre-sampling of the xfoosc func- 
tion. An MCMC final sampling is subsequently 
used to determine the 68% and 95% confidence- 
level contours for each constraint scenario stud- 
ied 4. 

For example, in [4] we showed that the ef- 
fect of dropping the f2cDM experimental data 
from the fit is not very important in constrain- 
ing the allowed regions in the ('tI]^/27"^o) a-nd 
(too, tan /3) planes. The reason for this can be 
understood by recalling that the WMAP strips 
in the CMSSM (toi/2,too) planes found for dif- 
ferent, but fixed, values of tan/3 move around as 
this and other CMSSM parameters are varied. In- 
deed, for fixed Aq, the strips can be shown to 
nearly foliate the (toi/2,'71o) plane [19,65]. Since 
tan P is only weakly constrained by the experi- 
mental data but gets correlated through the fit 
to the other parameters {mi/2, "tio, Aq), the ef- 
fect of the r^cDM constraint is to reduce the di- 
mensionality of the allowed parameter-space to 
a certain "hyper-sheet" which, when viewed by 
fixing tan/3 to a particular value (i.e. slicing the 
sheet along the tan/3-axis), reduces to the ob- 



*We note that for parameter space regions having low 
probability density, statistical fluctuations can appear in 
the form of an "archipelago of islands" near the 95% con- 
fidence levels. Such statistical fluctuations simply reflect 
the lower MCMC sampling density in regions of low prob- 
ability. 



served strips in the {mi/2,mQ) planes. However, 
since this sheet is generally embedded in the full 
parameter space hyper-volume and is not diago- 
nalized along some particular parameter axis, a 
large range of values for (toq, toi/2, ta'ii/3) re- 
main statistically probable when considering the 
global fit and, from a strict statistical considera- 
tion, there are no strips of preferred regions. 

When we apply here a similar analysis to the 
{g — 2)^ constraint, we find a very different pic- 
ture. We exhibited already in [4] the effect of 
relaxing this constraint by some fraction, show- 
ing that the preferred areas of the (to2/2: 1^0) and 
(toq, tan (3) planes changed substantially. Here we 
illustrate the effect of removing the {g — 2)^; con- 
straint entirely, which relaxes very considerably 
the upper limits on sparticle masses. However, 
the other observables still disfavour very large val- 
ues of Too and mi/2 by Ax^ ~ 2, as we discuss 
below. 

3. Summaries ofthe CMSSM and NUHMl 
Analyses 

The experimental constraints used in our anal- 
yses are listed in Table 1. The notations for 
the observables are standard, and were defined 
in [53,4]. Their values generally have only minor 
updates from the values quoted there, but one 
important comment concerns our implementation 
of the LEP constraint on Mji. The value quoted 
in the Table was derived within the SM, and is 
applicable to the CMSSM, in which the relevant 
Higgs couplings are very similar to those in the 
SM [66,67], so that the SM exclusion results can 
be used, supplemented with an additional theo- 
retical uncertainty whose implementation we now 
describe. 

We evaluate the x^{Mh) contribution within 
the CMSSM using the formula 

^ ^ (1.1 GeV)2 + (1.5 GeV)2 ' ^ 

with Mj^™" = 115.0 GeV for Mh < 115.0 GeV ^ 
Larger masses do not receive a x^iMh) contri- 

^We use 115.0 GcV so as to incorporate a conservative 
consideration of experimental systematic effects. 
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Observable 


Th. Source 


Ex. Source 


Constraint 


Add. Th. Unc. 


mt [GeV] 


[68,69] 


[70] 


173.1 ±1.3 






[68,69] 


[71] 


0.02758 ±0.00035 




Mz [GcV] 


[68,69] 


[71] 


91. 1875 ± 0.0021 




Ty fGcVl 

^ Zi [vJ"^ V J 


f68 691 


f711 


2.4952 ± 0.0023 


0.001 


crP , fnbl 

"had P^J 


f68 691 


f711 


41 540 ± n 037 




i?/ 


[68,69] 


[71] 

L ' J 


20.767 ± 0.025 






[68,69] 


[71] 


0.01714 ±0.00095 






[68,69] 


[71] 


0.1465 ± 0.0032 




i?b 


[68,69] 


[71] 


0.21629 ± 0.00066 


- 


i?c 


[68,69] 


[71] 


0.1721 ± 0.003 


- 


^fb(&) 


[68,69] 


[71] 


0.0992 ± 0.0016 


- 


Afb(c) 


[68,69] 


[71] 


0.0707 ± 0.0035 




^6 


[68,69] 


[71] 


0.923 ± 0.020 


- 


A, 


[68,69] 


[71] 


0.670 ± 0.027 


- 


A,(SLD) 


[68,69] 


[71] 


0.1513 ± 0.0021 


- 


sin^ 6*1, (Qfb) 


[68,69] 


[71] 


0.2324 ± 0.0012 




Mm/ [GeV] 


[68,69] 


[72,73] 


80.399 ± 0.025 


0.010 




[74-78] 


[79] 


1.117 ±0.076exp±0.082th(SM) 


0.050 


BR(B, ^ M+A^^) 


[80-83] 


[79] 


< 4.7 X 10-8 


0.02 X 10-8 




[82-84] 


[85-87] 


1.25±0.40[e,p+th] 






[80-83] 


[79] 


< 2.3 X 10-*^ 


0.01 X 10-9 




[88] 


[79,89] 


0.99 ±0.32 




BRk'^,./BR|^_,. 


[82,84] 


[90] 


1.008 ±0.014[e,p+th] 




T^Tjcxp /T3r)SM 
BR^^^^p/ BR^_^^^p 


[91] 


[92] 


< 4.5 




AM^^P/AAfSM 


[91] 


[93,94] 


0.97±0.01exp±0.27th(SM) 




(AA/-;/AA/™) 


[80-83] 


[79,93,94] 


1.00±0.01exp±0.13th(SM) 






[91] 


[93,94] 


1.08±0.14[exp+th] 






[95-98] 


[99-101] 


(30.2 ±8.8) X 10-1" 


2.0 X 10-1" 


M,. [GeV] 


[102-105] 


[106,107] 


> 114.4 (see text) 


1.5 




[108-110] 


[111] 


0.1099 ±0.0062 


0.012 



Table 1 



List of experimental constraints used in this work. The values and errors shown are the current best 
understanding of these constraints. The rightmost column displays additional theoretical uncertainties 
taken into account when implementing these constraints in the MSSM. We have furthermore taken into 
account the direct searches for SUSY particles at LEP [64]. 
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bution. The 1.5 GeV in the denominator corre- 
sponds to a convolution of the Ukelihood function 
with a Gaussian function. $i 5(2;). normahzed to 
unity and centered around M^, whose width is 
1.5 GeV, representing the theory uncertainty on 
Alh [102]. In this way, a theoretical uncertainty 
of up to 3 GeV is assigned for 95% of all Alh 
values corresponding to one CMSSM parameter 
point. The 1.1 GeV term in the denominator cor- 
responds to a parametrization of the CLg curve 
given in the final SM LEP Higgs result [106]. 

Within the NUHMl the situation is somewhat 
more involved, since, for instance, a strong sup- 
pression of the ZZh coupling can occur, inval- 
idating the SM exclusion bounds. In order to 
find a more reliable 95% C.L. exclusion limit for 
Mh in the case that the SM limit cannot be ap- 
plied, we use the following procedure. The main 
exclusion bound from LEP searches comes from 
the channel e+e^ ZH, H bb. The Higgs 
boson mass limit in this channel is given as a 
fimction of the ZZH coupling in [107]. A re- 
duction in the ZZh coupling in the NUHMl rel- 
ative to its SM value can be translated into a 
lower limit on the lightest NUHMl Higgs mass, 
M^™'*'°, shifted to lower values with respect to 
the SM limit of 114.4 GeV. (The actual num- 
ber is obtained using the code HiggsBounds [112] 
that incorporates the LEP (and Tevatron) lim- 
its on neutral Higgs boson searches.) For values 
of Mh < 86 GeV the reduction of the ZZh cou- 
plings required to evade the LEP bounds becomes 
very strong, and we add a brick-wall contribution 
to the function below this value (which has 
no influence on our results). Finally, eq. (2) is 
used with A/1'""* = a4'""*'° -I- 0.6 GeV to ensure 
a smooth transition to the SM case, where we 
use Mji'""' = 115.0 GeV to allow for experimen- 
tal systematics, as discussed above. This is a con- 
servative approach in the sense that the 1.1 GeV 
term used in eq. (2) can be regarded as a lower 
limit on the spread of the CLg curve in the vicin- 
ity of M;^™"'°. 

The numerical evaluation of the frequentist 
likelihood function using these constraints has 
been performed with the MasterCode [53,4], 
which includes the following theoretical codes. 



For the RGE running of the soft SUSY-breaking 
parameters, it uses SoftSUSY [113], which is com- 
bined consistently with the codes used for the 
various low-energy observables. At the elec- 
troweak scale we have included various codes: 
FeynHiggs [102-105] is used for the evaluation 
of the Higgs masses and (optionally) a^^^^ (see 
also [95-98])^. For flavour-related observables we 
use SuFla [82,83] as well as Superlso [119,120], 
and for the electroweak precision data we have in- 
cluded a code based on [68,69]. Finally, for dark- 
matter-related observables, MicrOMEGAs [108- 
110] and DarkSUSY [121,122] have been used. We 
made extensive use of the SUSY Les Houches Ac- 
cord [123,124] in the combination of the various 
codes within the MasterCode. 

It is well known from previous comparisons 
that different RGE codes for the running of the 
soft SUSY-breaking parameters give quite differ- 
ent results in parameter regions where higher- 
order corrections get very large [125]. This hap- 
pens in general for very large values of tan /9, but 
instabilities can also occur in, e.g., the evalua- 
tion of Ma in the CMSSM. In such a case the 
evaluation of the impacts of constraints that are 
affected by the heavy Higgs bosons can become 
unreliable. Motivated by these observations, we 
made two cuts on the pre-sampled CMSSM and 
NUHMl points in deriving the results presented 
below: we do not consider points with tan/3 > 60, 
and we have discarded parameter points where 
the difference between the running (DR) mass of 
the pseudoscalar Higgs, Ma{Q), and the physi- 
cal (on-shell) mass, Ma, gets unacceptably large. 
For the latter, we have applied the condition 
\y/Ml{Q) - Ma\/Ma > 0.6. Imposing these 
cuts has no effect on the best- fit point, nor on 
the 68% C.L. range of any parameter of the fit ^. 



^We recall that the experimental value appears to differ by 
over three standard deviations from the best SM calcula- 
tion based on low-energy e+e" data [95,99,100,114—117], 
but that the discrepancy is significantly reduced if r decay 
data are used to evaluate the SM prediction. We note that 
recently a new t based analysis has appeared [118], which 
yields a 1.9(7 deviation from the SM prediction. A new 
SM prediction based on radiative-return data from BaBar 
is in the offing. 

^However, discarding these points does reduce the 95% 
C.L. upper limit on nii^^ i^i the NUHMl by about 10%, 
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Motivated by {g — 2)^ and (to a lesser extent) 
BR(6 — > 57), we restrict our study to fi > 0. 

For the parameters of the best-fit CMSSM 
point we find tuq = 60 GeV, toi/2 — 310 GeV, 
Ao = 130 GeV, tan/3 = 11 and fi = 400 GeV, 
yielding the overaU xV^dof = 20.6/19 (36% 
probability) and nominally Mh = 114.2 GeV^. 
These values are very close to the ones reported 
in [4] . The corresponding parameters of the best- 
fit NUHMl point are mo = 150 GeV, mi/2 = 
270 GeV, Ao = -1300 GeV, tan/3 = 11 and 
ml^ = mf^^ = — 1.2 X 10^ GeV^ or, equivalently, 
^ = 1140 GeV, yielding = 18.4 (correspond- 
ing to a similar fit probability to the CMSSM) 
and Mji = 120.7 GeV. The similarity between 
the best-fit values of mp, mi/2 and tan/3 in the 
CMSSM and the NUHMl suggest that the model 
frameworks used are reasonably stable: if they 
had been very different, one might well have won- 
dered what would be the effect of introducing ad- 
ditional parameters, as in the NUHM2 [126,127] 
with two non-universality parameters in the Higgs 
sector^. 

These best-fit points are both in the coannihi- 
lation region of the {rriQ, 1111^2) plane, as can be 
seen in Fig. 1, which displays contours of the Ax^ 
function in the CMSSM (left) and the NUHMl 
(right). The C.L. contours extend to slightly 
larger values of mg in the CMSSM, while they 
extend to slightly larger values of mi/2 in the 
NUHMl, as was already shown in [4] for the 68% 
and 95% C.L. contours. However, the qualitative 
features of the Ax^ contours are quite similar in 
the two models, indicating that the preference for 
small niQ and mi/2 are quite stable and do not de- 
pend on details of the Higgs sector. We recall that 
it was found in [4] that the focus-point region was 
disfavoured at beyond the 95% C.L. in both the 
CMSSM and the NUHMl. We see in Fig. 1 that 
this region is disfavoured at the level Ax^ ~ 8 in 

from ~ 1000 GeV to ~ 900 GeV. This diflference may 
be regarded as a theoretical systematic uncertainty in the 
results. 

*This is acceptable, taking into account the theoretical 
uncertainty in the FeynHiggs calculation of A/j, [102] , see 
the discussion above. 

^Computationally, exploring adequately the NUHM2 pa- 
rameter space using the frequentist approach would be 
very expensive, but we hope to return to it in the future. 



the CMSSM and > 9 in the NUHMl. 

This feature is seen explicitly in the left and 
right panels of Fig 2, which display the likelihood 
functions for mo in the CMSSM and NUHMl, re- 
spectively. (We recall that the focus-point region 
would be found at niQ > 1500 GeV.) Looking 
first at the solid lines corresponding to the full 
global fit, we also see explicitly that low values of 
mo ~ 100 GeV are favoured in both cases, reflect- 
ing the fact that coannihilation points are gener- 
ally favoured. The favoured regions also have rel- 
atively low values of mi/2, as seen in Fig. 1. As 
we discuss in more detail later, the minimum in 
both cases is found at low tan /3 ~ 1 1 . 

The large values of Ax^ in the focus-point 
region are largely, but not entirely, due to the 
(g — 2)^ constraint, as can be seen in the dashed 
lines in Fig. 2, where this constraint has been re- 
moved. In the CMSSM case without (5-2)^, the 
global minimum at mo ~ 100 GeV is followed by 
a local maximum around mo ^ 1000 GeV with 
~ 3. This is in turn followed by a sec- 
ondary local minimum around mo ^ 2000 GeV 
with Ax^ ^ 2. The absolute minimum occurs 
in the coannihilation region, and the secondary 
minimum occurs in the focus-point region. The 
local maximum at intermediate mo reflects the 
fact that such values of mo are compatible with 
the dark matter constraint only at relatively large 
values of tan/3 and mi/2 that are disfavoured by 
other constraints. This is not the case in the 
NUHMl, where intermediate values of mo with 
relatively low values of tan /3 are compatible with 
the U^h^ constraint (thanks to the possible ap- 
pearance of direct-channel Higgs poles), as well 
as the other constraints. See Section 6 below for 
a more detailed discussion of the impact of drop- 
ping the {g — 2)^ constraint. 

We summarize in Table 2 the contributions 
to the global x^ likelihood function at the best- 
fit points in the CMSSM and NUHMl due to 
the most important observables as well as their 
total x^- We also list the contributions to x^ 
for the best fit we find in the focus-point (FP) 
region for the CMSSM (considered to be that 
with mo > 1000 GeV). This point has mo = 
2550 GeV, 7711/2 = 370 GeV, Ao = 1730 GeV 
and tan/3 = 51. It is apparent from Table 2 that 



irip [GeV/c^] 



m„ [GeV/c^] 



Figure 1. The Ax^ functions in the (mo,mi/2) planes for the CMSSM (left plot) and for the NUHMl 
(right plot). We see that the coannihilation regions at low mo and ?7ii/2 are favoured in both cases. 
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Figure 2. The likelihood functions for mo in the CMSSM (left plot) and in the NUHMl (right plot). The 
values are shown including (excluding) the {g ~ 2)^ constraint as the solid (dashed) curves. 



the focus-point region is disfavoured by {g — 2)^, 
but also by Mw, and that the contributions of 
the other observables fail to overcome this dis- 
advantage. Indeed, many of the other observ- 
ables favour independently the coannihilation re- 
gion, e.g., BR{Bu TVr), ^KSLD) and 



^"We note, however, that within the SM there is signifi- 
cant tension between the experimental value of (SLD) 



— though the difference here is relatively small. 
Only Mh and BR(6 57) and Afb(b)(LEP) 
favour the FP region, but not with high signif- 
icance. 



and Af[,(6)(LEP), and that this tension is not reduced sig- 
nificantly in the CMSSM or NUHMl, see also [69]. 
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Observable 


Best CMSSM fit 


Best NUHMl fit 


Best CMSSM FP fit 


(.9 - 2)^ 


0.44 


0.002 


8.4 




0.20 


0.41 


0.85 


Mw 


0.53 


0.08 


1.5 


Ae{SLD) 


2.84 


3.22 


3.56 


Afb(6)(LEP) 


7.61 


7.08 


6.74 


Ri 


0.96 


1.01 


1.05 


p,T>SUSY /rrSM 


1.16 


0.001 


0.95 


Mh 


0.17 








Xtot 


20.6 


18.5 


29.8 



Table 2 

The principal contributions to the global likelihood function from the experimental and phenomeno- 
logical constraints used in this work, as well as the total , for the best- fit points in the CMSSM and 
NUHMl, which both lie in the coannihilation region. For comparison, we also show the analogous num- 
bers for the best CMSSM fit we find in the focus-point (FP) region with niQ > 1000 GeV. Only those 
observables yielding the main contributions to the total are listed in the table. 



4. Likelihood Distributions for Sparticle 
Masses and Other Observables 

In our previous paper [4] we discussed, in ad- 
dition to the spectra at the best-fit points in 
the CMSSM and NUHMl, the regions of the 
(mo,TOi/2) planes preferred in these scenarios at 
the 68 and 95% C.L. Here we complement those 
discussions by providing directly the likelihood 
functions for certain sparticle masses, noting in 
particular the impacts of the most relevant con- 
straints^^. 

We start by discussing the likelihood functions 
for the mass of the neutralino LSP, m^o , in the 
CMSSM and NUHMl. The left panel of Fig. 3 
displays the likelihood function in the CMSSM. 
The solid line shows the result obtained when 
incorporating the LEP Higgs limit, while the 
dashed line corresponds to the case where the 
LEP Higgs constraint is removed. There is a 
sharp rise in the likelihood function at low val- 
ues of m^o, which is caused by the limits from 
the direct searches for SUSY particles, but re- 
ceives also contributions from BR(& — > 57) and 

^'^In each case, we show Ax^, the diflference between the 
total function and its value at the minimum for the 
relevant model. 



other constraints. This sharp rise in the likeli- 
hood function persists when the LEP Higgs con- 
straint is removed, but is shifted towards slightly 
lower values of m-^o in that case. The right panel 
of Fig. 3 shows the likelihood function for m^o 
in the NUHMl, again with and without the LEP 
Mh constraint imposed. Including the LEP Mh 
constraint we see that the optimal value of m^o 
is somewhat smaller than in the CMSSM case, 
reflecting the lower value of mi/2 at the corre- 
sponding best-fit point discussed in [4] Fi- 
nally, the dashed line in the right panel of Fig. 3 
displays the likelihood function in the NUHMl 
with the LEP Higgs constraint removed. Here we 
see very little difference from the result for the 
NUHMl with the LEP constraint imposed. This 
reflects the fact that in the NUHMl (unlike the 
CMSSM) the other constraints do not push Mh 
down to quite low values, a point made explicit 
in Fig. 4 below. 

The gradual rises in the likelihood functions at 
large m^o in both the CMSSM and the NUHMl 
arc dominated by the contribution of (g— 2)^, dis- 



^^We recall that, to a very good approximation, m^a ~ 
0.42 miy2 in most of the relevant regions of the CMSSM 
and NUHMl parameter spaces discussed here. 
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Figure 3. The likelihood functions for jti^o in the CMSSM (left) and in the NUHMl (right), both with 
(solid lines) and without (dashed lines) the LEP constraint on Mh- 



cussed already above, which disfavours large mi/2 
(and Too)- We comment later on the impacts on 
TOjjo and other observables if the {g — 2)^ con- 
straint is removed. 

In order to see explicitly the importance of the 
Mh constraint, we display in Fig. 4 the likeli- 
hood functions for Mh in the CMSSM (left) and 
the NUHMl (right), both with (solid lines) and 
without (dashed lines) the LEP constraint on Mh- 
Comparing first the two CMSSM results, we see 
that the other constraints would prefer a value 
of Mh somewhat below the SM Higgs limit from 
LEP [106] (this was already observed in [53]). 
The best fit value for Mh is still acceptable in 
that case, in particular in view of the theoretical 
uncertainties in the CMSSM evaluation of Mh, 
see the discussion above. However, in the case of 
the NUHMl, shown in the right plot of Fig. 4, the 
best-fit value of Mh indicated by the other con- 
straints is significantly higher than the SM LEP 
lower limit. As a consequence, incorporating the 
LEP constraint (sec above for details), shown as 
the solid line, does not alter significantly the best- 
fit value of Mh- As a corollary, the differences 
between the likelihood functions of the NUHMl 
for other masses and observables between the fits 



with and without the LEP Mh constraint are less 
significant than for the CMSSM. In the rest of 
this paper (except in Section 6.4) we show results 
with the LEP Mh constraint imposed. 

We discuss next the likelihood functions for 
various sparticle masses, which are summarized 
in Fig. 5. The results for the CMSSM spectrum 
are shown in the left plot, and for the NUHMl 
in the right plot. We start our discussion with 
the gluino mass, rUg. In both the CMSSM and 
the NUHMl, the best-fit points have relatively 
low values of rUg ~ 750 and ~ 600 GeV, respec- 
tively. These favoured values are well within the 
range even of the early operations of the LHC 
with reduced centre-of-mass energy and limited 
luminosity. However, the effect of the gradual in- 
crease in as TOi/2 increases, due essentially to 



{g — 2)^ as commented before, means that even 
quite large values of rUg ^ 2.5 TeV are allowed 
at the 3-cr (Ax^ = 9) level (not shown in Fig. 5). 
The LHC should be able to discover a gluino with 
rUg ~ 2.5 TeV with 100/fb of integrated luminos- 
ity at = 14 TeV [128,129], and the proposed 
SLHC luminosity upgrade to 1000/fb of inte- 
grated luminosity at ^/s = 14 TeV should permit 
the discovery of a gluino with rrig ~ 3 TeV [130]. 
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Figure 4. The likelihood functions for in the CMSSM (left) and in the NUHMl (right), both with 
(solid lines) and without (dashed lines) the LEP constraint on Mh- 



However, Fig. 5 does demonstrate that, whilst 
there are good prospects for discovering SUSY 
in early LHC running [4], this cannot be 'guaran- 
teed', even if one accepts the {g — 2)^ constraint. 

The central values of the masses of the su- 
persymmetric partners of the u, d, s, c, b quarks 
are slightly lighter than the gluino, as seen in 
Fig. 5. The difference between the gluino and the 
squark masses is sensitive primarily to mo. Since 
the preferred regions of the parameter space in 
both the CMSSM and the NUHMl are in the Xr 
slepton coannihilation region where mp < toi/2, 
mo makes only small contributions to the cen- 
tral values of the squark masses The SUSY 
partners of the left-handed components of the 
four lightest quarks, the ql, are predicted to 
be slightly heavier than the corresponding right- 
handed squarks, qr, as seen by comparing the 
mass ranges in Fig. 5. As in the case of the gluino, 
squark masses up to 2.5 TeV arc allowed at the 
3-cr level. Comparing the left and right panels, we 
see that the squarks arc predicted to be somewhat 
lighter in the NUHMl than in the CMSSM, but 
this difference is small compared with the widths 



However, this is not true in general, as we discuss in 
more detail later. 



of the corresponding likelihood functions. 

Turning now to the likelihood functions for the 
mass of the lighter stop, m^^, we find that it is 
shifted to values somewhat lower than for the 
other squark flavours. It can also be seen that 
the 2-cr range of its likelihood function differ from 
those of the gluino and the other squarks, reflect- 
ing the importance of scalar top mixing. We re- 
call that this depends strongly on the trilinear 
soft SUSY-breaking parameter At and the Higgs 
mixing parameter fi, as well as on the precise 
value of mt- As we discuss below, the favoured 
range of values of fj, is quite circumscribed in the 
CMSSM, whereas a larger variation in is possi- 
ble in the NUHMl. This has the effect of some- 
what broadening the likelihood function for 
in the NUHMl. 

In the case of the lighter stau fi, see its range 
in Fig. 5 and its likelihood ftmction in Fig. 6, the 
mass is very similar to that of the LSP Xi in the 
coannihilation region, but this is not the case in 
the rapid-annihilation H, A funnel region. The 
differences in the likelihood functions for the fi 
and the LSP Xij shown in Fig. 3, reflect the im- 
portance of this funnel region. In the case of the 
CMSSM (left panel of Fig. 6), the funnel region 
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Figure 5. Spectra in the CMSSM (left) and the NUHMl (right). The vertical solid lines indicate the 
best-fit values, the horizontal solid lines are the 68% C.L. ranges, and the horizontal dashed lines are the 
95% C.L. ranges for the indicated mass parameters. 



appears only at large values of tan (3 that are rela- 
tively disfavoured. This is why the shape of the fi 
likelihood function differs significantly from that 
of the Xi only at relatively large masses. In the 
case of the NUHMl (shown in the right panel of 
Fig. 6), rapid annihilation is possible also for low 
tan/?, leading to larger values of mo than in the 
CMSSM also for relatively small values of m^^ . 
The scalar taus as well as the other scalar lep- 



tons are expected to be relatively light, as can be 
seen in Fig. 5. They would partially be in the 
reach of the ILC(500) (i.e. with = 500 GeV) 
and at the 95% C.L. nearly all be in the reach of 
the ILC(IOOO) [131,132]. This also holds for the 
two lighter neutralinos and the light chargino (In 
the NUHMl, small parts of the the 95% C.L. re- 
gions for the masses of the heavier stau and the 
light chargino are above 500 GeV.) 
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Figure 6. The likelihood functions for nif-^ in the CMSSM (left panel) and in the NUHMl (right panel). 



The left plot of Fig. 7 displays the likelihood 
functions for fj, in the CMSSM (solid lines) and 
the NUHMl (dashed hues) In the CMSSM, 
the values of and Al^ are fixed in terms of 
the other model parameters by the electroweak 
boundary conditions. Consequently, the range of 
values for /i is quite small in the CMSSM, and 
the magnitude of fx turns out to be relatively 
small. In the NUHMl, the much larger range 
of /i reflects the greater freedom in the Higgs 
sector. Solving the electroweak vacuum condi- 
tions for models with non-universal Higgs masses 
broadens the /i distribution^^, with the implica- 
tions discussed above for the likelihood function 
for the ti. The right panel of Fig. 7 displays the 
likelihood functions for Ma (see also the range in 
Fig. 5). The likchhood function in the CMSSM 
is again somewhat narrower than in the NUHMl, 
reflecting the influence of the electroweak bound- 
ary conditions. The best-fit value in the CMSSM 
is significantly higher than in the NUHMl. Val- 
ues up to Ma ^ 500 GeV could be tested at the 
ILC(IOOO), i.e. the preferred regions of both mod- 
els could be probed. 

Fig. 8 displays the likelihood functions for 
tan/3. These are largely similar in the CMSSM 



^'^Wc recall that, motivated by (g — 2)^, and BR(6 S7), 
we study only fi > 0. 

^^Very large values of \fi\ > 1 TeV are disfavoured by the 
presence of deep charge- and colour-breaking minima [133, 
134], but this constraint is not applied here. 



and the NUHMl, with tan /? ^ 11 being favoured 
in both models. 

We turn now to the predictions for two other 
observables, namely BR(_Bs fj,~^^~) shown in 
Fig. 9 and the spin-independent Xi-proton scat- 
tering cross section a^^ shown in Fig. 10 
We see in the left panel of Fig. 9 that values 
of the BR{Bs fj,~^iJ.~) similar to that in the 
SM are favoured, particularly for the preferred 
lower values of tan/?. However, large deviations 
from the SM prediction (indicated by the vertical 
lines, which include the theoretical uncertainty) 
are still possible at the 3-cr level. The picture in 
the NUHMl is completely different, since the 
function is quite flat, with no significant penalty 
for substantial deviations from the SM prediction, 
and very large values of the branching ratio be- 
ing allowed at the 2-a level. The difference is 
largely due to the fact that smaller masses of 
the heavier Higgs bosons are permitted in the 
NUHMl. A large value of BR(S, ^ A'+M") 
would be a promising harbinger of SUSY at the 
LHC, and would favour a priori the NUHMl 
over the CMSSM. Assuming the SM value, i.e. 
BR{Bs -> M^/i") ~ 3.4 X 10"^, it has been es- 
timated [135] that LHCb could observe this pro- 
cess at the 5cr level within a few years of running. 



°The spin-independent Xi-proton and -neutron scatter- 
ing cross sections are very similar, and the spin-dependent 
scattering cross sections (not shown) are much further 
away from the prospective experimental sensitivity. 
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Figure 7. The likelihood functions for fi (left panel) and Ma (right panel) in the CMSSM (solid lines) 
and in the NUHMl (dashed lines). 
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Figure 8. The likelihood functions for tan (3 in the CMSSM (left panel) and in the NUHMl (right panel). 



This makes this process a very interesting probe 
of SUSY that could help to distinguish between 
different models. 

The value of (j^^ shown in Fig. 10 is calcu- 
lated assuming a tt-N scattering a term S^y = 
64 MeV: plausible values range between about 45 
and 80 MeV, and increases quite rapidly with 
Ejv [136-138]. We see in Fig. 10 that values of 
the X]'"P''oton cross section (7^^ ~ 10~^ pb are 
expected in the CMSSM, and that much larger 



values seem quite unlikely. On the other hand, in 
the NUHMl, though the best-fit value of the cross 
section is somewhat lower, a much larger range 
is possible Hence, detection of dark matter 
with a cross section much larger than ~ 10^® pb 
= lO"**^ cm^ would also be a good diagnostic 
for discriminating between the NUHMl and the 



^ 'No scaling of the cross-section was done here to account 
for regions where Q^h^ falls below the WMAP range, but 
such points pay a penalty. 
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Figure 9. The likelihood functions for the branching ratio BR(Bs /i^/i ) in the CMSSM (left panel) 
and in the NUHMl (right panel). The vertical lines indicate the SM value with its theoretical error. 



CMSSM. The present best upper limits on 
from the CDMS [139] and XenonlO [140] exper- 
iments arc at the ^ 10~^ pb level and the 
planned experiments should be sensitive down to 
below the - pb level [141,142]. 

5. Correlations between Sparticle Masses 
and with other Observables 

We now discuss in more detail some of the 
correlations between sparticle masses and observ- 
ables, starting with the LSP mass, rn^o, and the 
gluino mass, m^, shown in Fig. 11 . We ex- 
pect a very strong correlation, since the value of 
mi/2 largely controls both masses. However, in 
both cases there are radiative corrections that en- 
ter when making the transition from the SUSY- 
breaking parameters defined using the DR pre- 
scription to the on-shell masses, that depend on 
the other MSSM parameters. Moreover, the LSP 
is not a pure Bino, and the mixing with other 
neutralino states depends on the value of /i, in 
particular. Indeed, we see in Fig. 11a very strong 

Assuming a local LSP density of 0.3 GcV/cm'^, which is 
subject to astrophysical uncertainties. 

^^For one-dimcnsional scans, we continue to quote up 
to 9 units in Ax^, which corresponds to 3(T. For two- 
dimensional confidence level contour plots, we quote 1 — 
CL instead of the Ax^; the blue (red) lines in the plots 
correspond to 1 — CL = 32(5)%, and the white regions 
correspond to 1 — CL < 1%, or Ax^ > 9.21 units. 



m^o — rrig correlation in the CMSSM (left panel), 
which is not quite so strong for the NUHMl (right 
panel). Moreover, in the latter case we notice a 
small (grey) island of parameters where rrij^o is 
substantially lower than one would have expected 
for the corresponding value of m^. These few ex- 
amples have a Higgsino-like LSP, and have rela- 
tively small likelihoods. 

A corollary of the correlation between rrig and 
m^o seen in Fig. 11 is the relation between the 
mass scale of the heavy supersymmctric parti- 
cles [143] that might be discovered at the LHC 
with the threshold for producing the lighter spar- 
ticlcs that might be measured at a future linear 
e'^e~ collider. If one observes a gluino at a cer- 
tain mass scale (or establishes a lower limit on its 
mass), according to Fig. 11 one will have, within 
the CMSSM or the NUHMl, a lower bound on 
the threshold for pair-producing observable spar- 
ticles at a linear collider, y/s > 2m-^o The 
relevant m^o may be read directly off the verti- 
cal scale of Fig. 11. This is in general related to 
rrig by a simple, universal numerical factor, the 
only exception being the small island (which has 
a rather low likelihood) of models with unusually 
low m^o in the NUHMl, mentioned earlier and 
seen in the right panel of Fig. 11. 

In principle, the masses of the squark partners 

^''The lightest neutralino might then be visible in the 
channel e+e~ — » XiXi7 [144-146]. 
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Figure 10. The likelihood functions for the spin-independent x\-proton scattering cross section (in 
cn?) in the CMSSM (left panel) and in the NUHMl (right panel). 
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Figure 11. The correlation between the LSP mass, rn^o, and the gluino mass, nrig, in the CMSSM (left 
panel) and in the NUHMl (right panel). 



of the five lightest quarks depend on mo as well 
as mi/2. However, as seen in Fig. 12, they are 
also very highly correlated with nig, reflecting 
the fact that toq < mi/2 in the favoured regions 
of the CMSSM and the NUHMl, and also the 
fact that the sensitivities of m^^ ^ to toq are in- 
trinsically smaller than that to mi/2. That said, 
we see that the correlations of rUq^ with nig 
are slightly weakened in the CMSSM (left pan- 
els) at large rUg, reflecting the appearance of the 



rapid-annihilation funnel with relatively large mo 
at large mi/2 and tan/3. The greater width of 
the correlations in the NUHMl (right panels) at 
small rUg, compared to the CMSSM, reflects the 
possibility of greater toq due to the appearance 
of a rapid-annihilation funnel at smaller values of 
mi/2 and tan/3 than in the CMSSM. 

These effects are more visible in Fig. 13, where 
we plot the differences between the gluino and 
squark masses in the CMSSM (left plots) and in 
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Figure 12. The correlations between the gluino mass, nig, and the masses of the the left- and right-handed 
partners of the five light squark flavours, ruq^ ^ (upper and lower panels, respectively) are shown in the 
CMSSM (left panels) and in the NUHMl (right panels). 



the NUHMl (right plots). In the CMSSM, in the 
the cases of both the (upper left panel) and 
qR (lower left panel), we see that the squarks are 
always lighter than the gluino if rUg is itself light. 
However, if rrig > 1 TeV, although rUg > qL,R 
is still favoured, this is not necessarily the case, 
and rUg < qL,R becomes a possibility, because of 
the larger values of mo that occur in the rapid- 
annihilation funnel that appears as mi/2 and in- 
creases. In the case of the NUHMl (right panels), 
TTT-g < qL,R is & possibility also at low rug, thanks 
to the possible appearance of a rapid-annihilation 
funnel also at low mi/2. 

Fig. 14 displays the correlation between mg and 
TOj^ , which is somewhat weaker than the correla- 



tion between rrig and the other squark masses. 
This is because, in addition to sharing the de- 
pendence on mo with the other squarks, mi_^ is 
sensitive, as commented earlier, to the value of fi 
as well as mi/2 and mg. We recall further that 
the preferred range of /x is broader in the NUHMl 
than in the CMSSM, which explains why in this 
model the preferred range of mi_^ is broader for 
intermediate values of rUg. 

Fig. 15 displays the correlation between 
and rUg, which is generally proportional to the 
LSP mass, as discussed earlier. The rrifj - rUg 
correlation is strikingly different in the CMSSM 
(left panel) and the NUHMl (right panel). The 
tight correlation in the CMSSM reflects the fact 
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Figure 13. The differences between he gluino mass, nig, and the masses of the the left- and right-handed 
partners of the five light squark flavors, nriq^ ^ (upper and lower panels, respectively) and in the CMSSM 
(left panels) and in the NUHMl (right panels). 



that the favoured part of the parameter space is 
in the x?-fi coannihilation region, where the x?- 
fi mass difference is very smah. On the other 
hand, in the NUHMl, as already commented, 
there are favoured regions away from the coanni- 
hilation region, where rapid annihilation through 
direct-channel H, A poles keeps the relic density 
within the WMAP range. 

Fig. 16 demonstrates explicitly the big con- 
trast between the behaviours of the fi - Xi 
mass difference in the CMSSM (left panel) and 
the NUHMl (right panel). We see that in 
the CMSSM small mass differences are always 
favoured, and are mandatory for LSP masses 
^ 200 GcV, whereas larger mass differences are 



possible for LSP masses ^ 200 GeV, as the 
rapid-annihilation funnel opens up. However, in 
the NUHMl large mass differences are possible 
for all LSP masses, particularly for LSP masses 
^ 200 GeV. This means that, whereas in the 
CMSSM the 'visible' fi pair-production threshold 
at the ILC may be only slightly higher than the 
'invisible' x\ pair-production threshold, it may 
be considerably higher in the NUHMl, namely 
mf, < 400 GeV at the 95% C.L.. This is a po- 
tentially crucial signature for distinguishing the 
NUHMl from the CMSSM. 

Fig. 17 displays the favoured regions in the 
[Ma, tan /3) planes for the CMSSM and NUHMl. 
We see that they are broadly similar, with 
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Figure 14. The correlation between TOj^ 
NUHMl (right panel). 



and the gluino mass, rrig, in the CMSSM (left panel) and in the 




Figure 15. The correlation between rrif-^ and the gluino mass, nig, in the CMSSM (left panel) and in the 
NUHMl (right panel). 



little correlation between the two parameters. 
(M^,tan/3) planes in certain benchmark scenar- 
ios have often been used in the past to analyze 
the prospects for discovering heavy Higgs bosons 
at the LHC [147,148]. Most of these analyses 
have been done in the context of scenarios that 
do not take the relic-density constraint into ac- 
count, for exceptions sec [34.35]. The Higgs dis- 



covery contours determined in the various bench- 
mark scenarios cannot directly be applied to the 
(M^,tan/3) planes in Fig. 17 displaying our fit 
results for the CMSSM and the NUHMl. In or- 
der to assess the prospects for discovering heavy 
Higgs bosons at the LHC in this context, we fol- 
low the analysis in [149], which assumed 30 or 
60 fb"^ collected with the CMS detector. For 
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evaluating the Higgs-sector observables including 
higher-order corrections we use the soft SUSY- 
breaking parameters of the best-fit points in the 
CMSSM and the NUHMl, respectively. We show 
in Fig. 17 the 5-cr discovery contours for the three 
decay channels H,A^ t^t~ jets (solid lines), 
jet -I- ji (dashed lines) and jet + c (dotted lines). 
The parameter regions above and to the left of 
the curves are within reach of the LHC with 
about 30 fb~^ of integrated luminosity. We see 
that most of the highest-CL regions lie beyond 
this reach, particularly in the CMSSM. At the 
ILC(IOOO) masses up to Ma < 500 GeV can be 
probed. Within the CMSSM this includes the 
best-fit point, and within the NUHMl nearly the 
whole 68% C.L. area can be covered. 

Wc display in Fig. 18 the correlations between 
Ma and mi/2 in the CMSSM and in the NUHMl. 
In the former case, the electrowcak boundary con- 
ditions fix Ma, and the effect is to force Ma > 
2m-^o. However, Ma becomes essentially a free 
parameter in the NUHMl, and values smaller 
than m^o become possible also. On the other 
hand, there is a narrow strip where Ma ^ 2m-^o 
which is disfavoured because there rapid direct- 
channel annihilation suppresses the relic density 
below the range preferred by astrophysics and 
cosmology. The points with Ma < 2TOj^o are a 
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in the 



qualitatively new possibility opened up within the 
NUHMl as compared to the CMSSM, and extend 
to relatively large values of 7711/2. 

Fig. 19 displays the correlation between tan /? 
and the BR(i?s — > /i^/i"). As seen previously, in 
the CMSSM the preferred values of the branch- 
ing ratio arc very close to the value in the SM, 
though somewhat larger values may occur at large 
tan /3, which however have a lower likelihood. The 
situation is completely different in the NUHMl, 
where much larger values of the branching ra- 
tio for BR{Bs — > fX^ n~) are possible, even if 
tan/3 ~ 10. This increase reflects the possibil- 
ity that AIa may be considerably smaller than 
in the CMSSM. The upper right corner of the 
NUHMl plot, i.e., simultaneous large tan/? and 
large BR(i3s /i+/i"), is disfavoured because it 
would give rise to values of BR(fe S7) that are 
too small. 

Fig. 20 displays the preferred range of the spin- 
independent DM scattering cross section a^^ (cal- 
culated assuming a tt-N scattering a term Sat = 
64 MeV) as a function of m^a. In the case of 
the CMSSM, we see that the expected range of 
(jp^ lies mainly between the present experimental 
upper limits (solid lines) [139,140], which start to 
touch the preferred region, and the projected sen- 
sitivity of the SuperCDMS experiment (dashed 



Figure 16. The correlation between the fi - Xi raass difference as a function of the LSP mass, 
CMSSM (left panel) and in the NUHMl (right panel). 
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Figure 17. The correlations between Ma and tan/3 in the CMSSM (left panel) and in the NUHMl (right 
panel). Also shown are the 5-a discovery contours for observing the heavy MSSM Higgs bosons H,A in 
the three decay channels H,A^ t~^t~ — > jets (solid line), jet + /i (dashed line), jet + e (dotted line) at the 
LHC. The discovery contours have been obtained using an analysis that assumed 30 or 60 fb^^ collected 
with the CMS detector [129,149]. 
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Figure 18. The correlations between Ma and mi/2 in the CMSSM (left panel) and in the NUHMl (right 
panel). 



line) [142], which should cover the preferred re- 
gion. As noted earher, these experimental con- 
straints were not applied in our analysis. The 
uncertainty in S^r and the astrophysical uncer- 



tainties in the local dark matter density (which 
are difficult to quantify), preclude including the 
value of (Tp' in the likelihood analysis presented 
here. This region is in good agreement for neu- 
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Figure 19. The correlation between branching ratio for BRfBg —f iX^ ^ ) and tan/3 in the CMSSM (left 
panel) and in the NUHMl (right panel). 



tralino masses between 100-300 GeV with that 
found in [63], where a recent scan (without hkc- 
Uhood information) was performed. 

As aheady commented, the range in the 
NUHMl is larger than in the CMSSM. We see 
in Fig. 20 that the larger cross-section values oc- 
cur, as expected, for small nrij^o^, in particular in 
the small island of Higgsino-like DM that appears 
close to the 95% C.L. for m^o < 100 GeV. If 
100 GeV < TO^o < 200 GeV, the allowed range 
of the cross section is larger than in the CMSSM 
because of the wider range of possible values of 
Ma as found in [63] for this neutralino mass 
range, and the present experimental sensitivity 
is already below the values of cr|^ found for some 
favoured NUHMl parameter values assuming the 
nominal values of Sjv and the local LSP den- 
sity. The smallest values of the cross section oc- 
cur when m^o > 200 GeV, in models close to the 
95% C.L. limh for the NUHMl, which have Ma 
larger than in the CMSSM. In general, we see 
that whereas the favoured values of 0-^^ are close 
to the present experimental upper limits [155] in 
both the CMSSM and the NUHMl, there is a 
greater possibility in the NUHMl that the cross 
section may lie beyond the projected sensitivity 
of SuperCDMS [142]. 



6. Dropping Constraints 

6.1. Dropping the (g — 2)^ Constraint 

We have stressed above that the results in the 
previous section are strongly dependent on the 
implementation of the (g — 2)^ constraint. In par- 
ticular, we have displayed in Fig. 2 above the like- 
lihood functions for mo in the CMSSM (left) and 
the NUHMl (right) both with the {g - 2)^ con- 
straint imposed (solid) and without it (dashed). 
We now discuss in more detail the effect of drop- 
ping the {g — 2)^ constraint completely, calcu- 
lating a new Xkiosc with no contribution from 
(ff-2)^. 

The function obtained for mp in the 
CMSSM without the {g — 2)^ constraint, shown 
in the left panel of Fig. 2, is much flatter than 
the corresponding function obtained with the 
{g — 2)^ constraint. Nevertheless, we see non- 
trivial features in the function. One is that the 
location of the CMSSM global minimum is very 
similar to the case with the (5 — 2)^ constraint ap- 
plied. We recall that the rise in the function 
at small mo is determined essentially by the Mh 
and BR(& — > 57) constraints, with [g — 2)^j not 
playing a role. However, it is perhaps surprising 
that the other constraints cause x^ to rise until 
mo ~ 1000 GeV, where Ax^ ^ 3. However, we 
see in Table 2 that, in addition to (5 — 2)^, sev- 
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Figure 20. T/ie correlation between the spin-independent DM scattering cross section a^^ (calculated 
assuming a n-N scattering a term E^r = 64 MeV) and rn^o in the CMSSM (left panel) and in the 
NUHMl (right panel). The solid lines [155] are the present experimental upper limits from CDMS [139] 
and XENON10[140], and the dashed line [155] indicates the projected sensitivity of the SuperCDMS 
experiment [14-2]. 



eral other constraints favour the best-fit CMSSM 
point over points with mo > 1000 GeV, including 
BR(i3„ —f Ti'r), Mw, Ai and Ri. Continuing to 
larger mo in the left panel of Fig. 2, we see that 
decreases again slightly, but that still Ax^ ^ 2. 

Similar features are seen in the function ob- 
tained for mo in the NUHMl without the (g — 2)^ 
constraint, shown in the right panel of Fig. 2. 
Again, the value of mo at the best-fit NUHMl 
point is very similar, whether (17 — 2)^ is included 
or not, and again Ax^ ^ 2 at large mo. However, 
there is no intermediate hump at mo ~ 1000 GeV 
analogous to that in the CMSSM, reflecting the 
greater freedom in the NUHMl to adjust param- 
eters so as to obtain a lower value of x^- 

A corollary of the observations in the previous 
paragraphs is that, at some level, the other con- 
straints favour a non-zero supersymmetric contri- 
bution to {g—2)fj_. This is indeed visible in Fig. 21, 
where we see the predicted values of the contribu- 
tions of supersymmetric particles to (17—2)^ in the 
CMSSM (left) and the NUHMl (right). We show 
the x^ functions only for positive contributions to 
(g — 2)^, since our points were all chosen to have 
/i > 0. Nevertheless, the fact that the minima of 
the x^ distributions are for A((g — 2)^) ^ is 
non-trivial, because it reflects the above observa- 



tion that large values of the sparticle masses are 
disfavoured, and the order of magnitude predic- 
tion for A{{g — 2)^) agrees with estimates based 
on low-energy e'^e^ data. 

6.2. Dropping the BR(6 — »■ sj) Constraint 

We display in Fig. 22 the effects on the CMSSM 
and NUHMl fits (left and right panels, respec- 
tively) of omitting the BR(& 57) constraint 
from the global fit, as obtained by calculating 
a new X^oog^ with no contribution from BR(fe 
s"f). In both models, we see that the predictions 
for BR(6 — > 57) based on the other constraints 
(solid lines) are not very precise. The best-fit val- 
ues for BR(6 — *■ 57) are in both models quite close 
to the SM and hence the experimental value, but 
the CMSSM permits much smaller values, and 
both larger and smaller values are allowed in the 
NUHMl with relatively small increases in x^- 

The converse statement is that applying the 
BR(6 — !■ 57) constraint does not impose a large 
price on the global minimum. This is appar- 
ent from Table 2, where we saw that BR(6 57) 
contributes about Ax^ ~ 1 to the total in 
the CMSSM and yields a negligible contribution 
in the NUHMl. There is no tension between 
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Figure 21. The functions for the supersymmetric contributions to {g — 2)^ in the CMSSM (left) and 
the NUHMl (right), as calculated using the other constraints except [g — 2)^ itself (solid line), and with 
all constraints included (dashed line). 




Figure 22. The functions for the ratio of the MSSM prediction over the SM prediction to BR(b sj), 
R{b S7) = BRf^^^^/BR^'^,^, m the CMSSM (left) and the NUHMl (right), as calculated using the 
other constraints except BR(b s"f) itself (solid line), and with all constraints included (dashed line). 



BR(6 sj) and the other constraints. 
6.3. Dropping the fl-^^h^ Constraint 

One of the most exciting predictions of the 
CMSSM and the NUHMl is the existence of a 
cold dark matter candidate in the form of the 



LSP, which we assume here to be the lightest 
neutrahno [150]. It is natural to take the next 
step, and ask whether these models predict a 
relic LSP density that is close to the experimen- 
tal value of the cold dark matter density. This 
density is determined with an accuracy of a few 
percent, and an comparable accuracy in the pre- 
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diction based on the other available experimental 
constraints will be difhcult. This will improve if 
(when) the LHC discovers SUSY and its param- 
eters are measured more accurately at an e'^e~ 
linear collider [65,151,152,57]. 

Nevertheless, calculating a new xf^og;, with no 
contribution from fl^h'^, it is interesting to see in 
Fig. 23 that both the CMSSM (left) and NUHMl 
(right) favour ranges of fl^^h^ values that include 
the measured values of the cold dark matter den- 
sity. In the case of the CMSSM, the prediction 
for n^/i^ (solid line) is within an order of magni- 
tude above and below the measured value at the 
level Ax^ < 4. This is also true in the NUHMl 
above the measured value, but the relic LSP den- 
sity could be two or more orders of magnitude 
below the measured value with Ax^ < 1- This 
is because there is a possibility that the relic 
density may be suppressed by rapid annihilation 
through direct-channel Higgs poles in the region 
of relatively low mi/2 and tan/3 in the NUHMl 
that is favoured by the other constraints, notably 
(5-2)^. 

The converse statement is that applying the 
^x^^ constraint also does not impose a large 
price on the global minimum. In fact, Ax^ ^ 1 
in both the CMSSM and the NUHMl, and there- 
fore has not been listed in Table 2. As in the case 
of BR{b s-f), there is no tension between ^^^^ 
and the other constraints. 

6.4. Dropping the Mh Constraint 

We have already commented on the effect on the 
likelihood function for m^o of dropping the LEP 
Mh constraint, see Fig. 3, and on the prediction 
for Mfi itself, see Fig. 4. We now discuss in more 
detail the likelihood functions for Mh within the 
CMSSM and NUHMl frameworks obtained when 
dropping the contribution to from the direct 
Higgs searches at LEP, shown in the left and right 
panels of Fig. 24, respectively. The left plot up- 
dates that for the CMSSM given in [53]. 

It is well known that the central value of the 
Higgs mass in a SM fit to the precision elec- 
troweak data lies below 100 GeV [153,73], but the 
theoretical (blue band) and experimental uncer- 
tainties in the SM fit are such that they are still 



compatible at the 95% C.L. with the direct lower 
limit of 114.4 GeV [106] derived from searches at 
LEP. In the case of the CMSSM and NUHMl, 
one may predict Mh on the basis of the underly- 
ing model parameters, with a l-cr uncertainty of 
1.5 GeV [102], shown as a red band in Fig. 24. 
Also shown in Fig. 24 are the LEP exclusion on 
a SM Higgs (yellow shading) and the ranges that 
are theoretically inaccessible in the supersymmet- 
ric models studied (beige shading) The LEP 
exclusion is directly applicable to the CMSSM, 
since the h couplings are essentially indistinguish- 
able from those of the SM Higgs boson [66,67], 
but this is not necessarily the case in the NUHMl, 
as discussed earlier in this paper. 

In the case of the CMSSM, we see in the left 
panel of Fig. 24 that the minimum of the func- 
tion occurs below the LEP exclusion limit. While 
the tension between the x^ function for Mh aris- 
ing from the CMSSM fit and the LEP exclusion 
limit has slightly increased compared to the ear- 
lier analysis performed in [53] , the fit result is still 
compatible at the 95% C.L. with the search limit, 
similarly to the SM case. As we found in the anal- 
ysis performed above, a global fit including the 
LEP constraint has acceptable x^- In the case of 
the NUHMl, shown in the right panel of Fig. 24, 
we see that the minimum of the x^ function oc- 
curs above the LEP lower limit on the mass of a 
SM Higgs. Thus, within the NUHMl the combi- 
nation of all other experimental constraints nat- 
urally evades the LEP Higgs constraints, and no 
tension between Mh and the experimental bounds 
exists. 

7. Conclusions 

We have presented in this paper detailed results 
from global fits to available experimental and cos- 
mological data within the CMSSM and NUHMl, 
using a frequentist approach. As already reported 
in [4], we find relatively small values of the key 
input SUSY-breaking parameters mi/2 and mo 
in both models. Moreover, the values for these 
parameters are quite similar in the two models, 

^^It is apparent that the current Tevatron exclusion [154] 
of a range between 160 and 170 GeV does not impact 
supersymmetric models. 
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Figure 23. The functions for the supersymmetric contributions to il^h^ in the CMSSM (left) and 
the NUHMl (right), as calculated using the other constraints except fl^h^ itself (solid line), and with all 
constraints included (dashed line). 
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Figure 24. The functions for Mh in the CMSSM (left) and the NUHMl (right), including the theoretical 
uncertainties (red bands). Also shown is the mass range excluded for a SM-like Higgs boson (yellow 
shading), and the ranges theoretically inaccessible in the supersymmetric models studied. 



indicating that the predictions are relatively ro- 
bust and do not depend strongly on the details of 
the Higgs sector. 

We have presented details of the likelihood 
functions for individual sparticle masses and the 
correlations between them. As noted in [4], 
the particle spectra are similar in the two mod- 



els, the most prominent differences being in the 
masses of the heavier Higgs bosons A, H and H^, 
which are lighter in the NUHMl, and the heav- 
ier neutralinos and chargino, which are lighter in 
the CMSSM. These differences reflect the greater 
freedom in varying the parameters of the Higgs 
sector in the NUHMl. The favoured values of 
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the particle masses in both models are such that 
there arc good prospects for detecting supersym- 
mctric particles even in the early phase of the 
LHC running with reduced centre-of-mass energy 
and limited luminosity and for observing super- 
symmetric particles and possibly the whole Higgs 
boson spectrum at a 1 TeV e+e" collider (the lat- 
ter refers in particular to the case of the NUHMl). 

We find striking correlations between the dif- 
ferent sparticle masses in both the CMSSM and 
the NUHMl. This reflects the fact that the domi- 
nant contributions to most of the sfermion masses 
are those due to mi/2, implying that most spar- 
ticle masses are tightly correlated with those of 
the gluino and the LSP. These correlations imply 
that, if the gluino is discovered at the LHC and its 
mass determined by a combination of kinematic 
and cross-section measurements, the predictions 
for the other sparticle masses within the CMSSM 
and the NUHMl could be refined considerably. 
In particular, the masses for colour-neutral sparti- 
cles such as the neutralino LSP and sleptons could 
be estimated more accurately, and hence also the 
energies of the corresponding thresholds in e+e^ 
annihilation within these models. For some of the 
correlations, most notably the difference between 
the LSP mass and the mass of the lighter stau, 
the pattern of the fit results in the NUHMl dras- 
tically differs from the one in the CMSSM. Mass 
correlations of observed supersymmetric particles 
could therefore provide very valuable information 
for distinguishing between different models. 

In addition to the sparticle masses, there are 
several other observables that could serve to con- 
strain (or provide evidence for) the CMSSM or 
the NUHMl. One observable that could discrim- 
inate between the CMSSM and the NUHMl, and 
might lead to an early discovery at the LHCb ex- 
periment, is BR{Bs In the CMSSM 
the rate for BR(i?s fi^^^) obtained from 
the fit is expected to be close to the SM value, 
whereas the value may be considerably larger in 
the NUHMl without reducing the goodness of the 
fit. 

A very exciting measurement would be that of 
the direct scattering of astrophysical cold dark 
matter particles. We find in both the CMSSM 
and the NUHMl that the favoured rate for spin- 



independent dark matter scattering lies quite 
close to the present experimental upper limit, 
though with larger uncertainties in the NUHMl. 
In view of the prospective improvements in the 
sensitivities of direct dark matter search exper- 
iments in the near future, they may be able to 
actually find the first indication of a supersym- 
metric particle before the LHC, though a combi- 
nation of astrophysical and collider measurements 
would be needed to pin down its SUSY nature. 

We have emphasized throughout this paper the 
sensitivity of our conclusions to the imposition of 
the {g — 2)^ constraint. This plays the dominant 
role in disfavouring large values of mi/2 and niQ 
and hence, in particular, the focus-point region 
of the CMSSM. In particular, BR(& ^ 57) plays 
no role in disfavouring the focus-point region. In- 
triguingly, however, some other observables seem 
slightly to prefer independently the coannihila- 
tion region, such as Mw and JiK{Bu tv-t). 
The net result is that the focus-point region is 
disfavoured by Ax^ ^ 2, even if the {g — 2)^ con- 
straint is dropped. Conversely, the other data 
provide a hint that the supersymmetric contribu- 
tion to {g—2)^ might be of comparable magnitude 
to the range required to reconcile the experimen- 
tal measurement of (g — 2)^ with the SM calcu- 
lation. 

We have also explored the effect of dropping 
from the global fit the experimental measurement 
of BR(& 37), and have shown that there is 
no conflict between this observable and the other 
constraints. We have shown as well that if ^l^h^ is 
dropped from the global fit, the other constraints 
favour — quite remarkably — a range within an 
order of magnitude of the astrophysical cold dark 
matter density, particularly within the CMSSM. 
These studies reveal no latent tensions between 
the data and either the CMSSM or NUHMl fit. 
Finally, we have discussed the impact of drop- 
ping the LEP Higgs constraint from the global 
fits. While in the CMSSM there is a slight ten- 
sion between the fit result and the direct search 
limit, similarly to the SM case, the NUHMl ac- 
tually favours a value for significantly above 
the LEP limit. The discovery at the LHC of a 
Higgs boson weighing more than 120 GeV would 
favour the NUHMl over the CMSSM. 
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Indirect constraints on supersymmetric model 
parameters are fine in their own way, and it is 
encouraging that there are no significant tensions 
in either the CMSSM or NUHMl fits. However, 
we hope that soon it will be possible to include 
in such fits some experimental measurements of 
physics beyond the SM. The fit results seem to 
indicate that there are good prospects for sparti- 
cle and Higgs-boson production at the LHC, but 
that there may also be good chances at a similar 
time scale to obtain evidence for cold dark matter 
scattering or for a discrepancy with the SM pre- 
diction for some other observable besides (g — 2)^, 
such as BR(i?s fi^ n~) or BR(i?„ tVt)- 
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Note added 

While we were finalizing this paper, the anal- 
ysis [57] has appeared. This uses a previous ver- 
sion of the MasterCode [4,53] to fit available data 
within the CMSSM and also adopts a frequentist 
Markov Chain Monte Carlo approach: the results 
of the analysis are very similar to ours. Our paper 
also compares current CMSSM and NUHMl fits, 
whereas [57] discusses the perspectives for fits us- 
ing future data. 
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